Gluonic phase versus LOFF phase in two-flavor quark matter 
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We study the gluonic phase in a two-flavor color superconductor as a function of the ratio of 
the gap over the chemical potential mismatch, A/S/j,. We find that the gluonic phase resolves the 
chromomagnetic instability encountered in a two-flavor color superconductor for A/ 5 fx < y2. We 
also calculate approximately the free energies of the gluonic phase and the single plane-wave LOFF 
phase and show that the former is favored over the latter for a wide range of coupling strengths. 

PACS numbers: 12.38.-t, ll.30.Qc, 26.60.+C 

It is widely accepted that sufficiently cold and dense 
quark matter is a color superconductor The most 
likely and, probably, the only place where color supercon- 
ductivity can exist in the universe is the interior of com- 
pact stars. Thus, studies of phases of quark matter under 
conditions realized in the bulk of compact stars (i.e., color 
and electric charge neutrality, and /3-equilibrium) have 
recently attracted a great deal of interest. The density 
regime of relevance for compact stars is up to a few times 
the normal nuclear density po — 0.16 fm _ . In this "mod- 
erate" density regime, the studies of QCD-motivated ef- 
fective theories are most useful and have revealed a rich 
phase structure @, H, 0| • 

One of the most striking features of neutral and j3- 
cquilibrated color-superconducting phases is unconven- 
tional cross-flavor Cooper pairing of quarks with the pos- 
sibility of gapless superconductivity, e.g., in the form of 
the gapless 2SC (g2SC) phase || or the gapless color- 
flavor-locked (gCFL) phase @. It was, however, quickly 
realized that the 2SC/g2SC phases suffer from a chro- 
momagnetic instability, indicated by imaginary Meissner 
screening masses of some gluons Q. In the 2SC phase, 
these instabilities occur when the ratio of the gap over 
the mismatch of the chemical potential, A/<5/i, decreases 
below a value y/2. Similar instabilities were found also 
in the gCFL phase |[. 

Resolving the chromomagnetic instability and clarify- 
ing the nature of the true ground state of dense quark 
matter are the most pressing tasks in the study of color 
superconductors. It was proposed that the chromomag- 
netic instability in two-flavor quark matter can be re- 
moved by the formation of a single plane-wave LOFF 
state uBl [ill [l2| (first studied by Larkin and Ovchin- 
nikov [LI, and Fulde and Ferrell [14[ in the context of 
solid state physics, and by Alford, Bowers, and Rajagopal 



15] for cold, dense quark matter), or a gluonic phase 
with vector condensation in the ground state 11611. (For 
a recent discussion of this issue, see also Refs. [13, [lj.) 
Alternatives include a mixed phase [l!| and, in the case 
of three-flavor quark matter, also phases with sponta- 
neously induced meson supercurrents [2(|. While the 
neutral LOFF state is free from the chromomagnetic in- 
stability in the weak-coupling regime [To| . this is, in fact, 
not the case for somewhat larger values of the coupling 
[2l| . At the same time, the gluonic phase can resolve 
the instability there. So far, however, the gluonic phase 
phase has been studied only around the critical point 



The gluonic phase and the LOFF phases are currently 
viewed as the most likely candidates for resolving the 
chromomagnetic instability and, thus, for the true ground 
state of two-flavor color-superconducting quark matter. 
(Here we exclude the possibility of phase separation [n| 
which may have limitations of its own and deserves a 
separate in-depth study.) In order to see which of the 
two proposed phases is actually preferred, one first has 
to extend the analysis of Ref. . 16j to a computation of the 
free energy away from the critical point A/Sfj, — y2, and 
then compare the results to the free energy of the single- 
plane LOFF state. 1 This is done in the present work. We 
qualitatively confirm the results of Ref. [l8| and extend 
them by (approximately) including the neutrality condi- 
tion and explicitly comparing the free energy of the glu- 
onic phase to that of the single plane- wave LOFF phase 
as a function of the coupling strength. 



In order to study various phases of two-flavor quark 
matter, we use a gauged Nambu-Jona-Lasinio (NJL) 
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Note that while the comparison with a multi-plane wave LOFF 
states would be more desirable, the corresponding free energy 
cannot be easily estimated within a microscopic approach. For 
current state-of-the-art calculations using an effective theory see 
Ref. [51. 
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model with massless up and down quarks: 



(1) 



where the quark field ip carries flavor = l,...Nf 
with Nf = 2) and color (a,0 = 1, . . . N c with N c = 3) 
indices, C is the charge conjugation matrix; (e) lk — e lk 
and (e b ) a P = e ba P are the antisymmetric tensors in flavor 
and color spaces, respectively. The covariant derivative 
and the field strength tensor are defined as 
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pa 

jjLV 



= d„A a u -d v Al + gf 



abc 4 6 



(2a) 
(2b) 



To evaluate loop diagrams we use a three-momentum cut- 
off A. Hence, the model has two phenomenological model 
parameters, the cutoff A and the diquark coupling Go- 
We use A = 653.3 MeV throughout this paper, but we 
consider Go as a free parameter. Henceforth, in order 
to specify the diquark coupling we use Ao which is the 
value of the 2SC gap at Sfi = (see below). 

In /3-equilibrated neutral 2SC/g2SC matter, the ele- 
ments of the diagonal matrix of quark chemical potentials 
/} are given by 



Hub 



M«9 = P> - <fy*> 

Hdg = M + 

Jl-Sfi- /i 8 , 
ji + Sfj,- /i 8 , 



with 



M = M- T + T , *M= T - 



(3a) 
(3b) 
(3c) 
(3d) 



(4) 



In a gauge theory, the self-consistent solution of the 
Yang-Mills equations requires background gauge fields 
[23| . These can be viewed as electric- and color-chemical 
potentials which ensure electric and color-charge neutral- 
ity of the system. Note that a generalization of this holds 
true even in the case of inhomogeneous phases. Then, of 
course, the corresponding fields would not be constant in 
space. Instead, they would have a constant central value 
contribution and, on top of it, a coordinate-dependent 
modulation describing color-electric fields induced by the 
inhomogeneities. 2 The constant contribution would take 
care of the global neutrality, while the modulation de- 
scribes the local field needed to prevent the local flow of 
currents. 



On the other hand, in NJL-type models without dy- 
namic gauge fields, one has to ensure electric and color- 
charge neutrality by introducing appropriate chemical 
potentials by hand [24|. In the case of the 2SC/g2SC 
phases, we only require an electron chemical potential 
/i e , and a color-chemical potential /is which ensures that 
the color-charge density ns is zero. In principle, in other 
phases like the gluonic phase one has to check that no 
other color-charge density is non-vanishing and necessi- 
tates the introduction of a respective color-chemical po- 
tential. Indeed, the gluonic phase introduced in Ref. 
[l6| requires a non- vanishing temporal component of the 
gluon field of the third color, (Aq). In our gauged 
NJL model, this is equivalent to a non-vanishing color- 
chemical potential fi3 besides /is ■ In this first exploratory 
study, however, we use the fact that both fi^ and /is are 
known to be numerically small and we simply neglect 
them. 

In Nambu-Gor'kov space, the inverse full quark prop- 
agator S -1 ^) is written as 



(So")- 1 



(5) 



with 



(S+y 1 = 7 % + {fi- <5/ir 3 ) 7 ° + <n»AlT*, (6a) 
(S^)- 1 = 7% - (/2 - */iT 3 )7° - grA-T aT ,(6h) 



and 



$ = -iee 75 A 



(7) 



Here t 3 = diag(l, —1) is a matrix in flavor space. Follow- 
ing the usual convention, we choose the diquark conden- 
sate to point in the third (blue) direction in color space. 

In the one-loop approximation, the free energy of two- 
flavor quark matter at T = is given by 



Va 



4G 



D 



d 4 p 
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lnDetS"- 1 ^), 



(8) 



where "Det" stands for the determinant in Dirac, flavor, 
color, and Nambu-Gor'kov space. Unlike the free energy 
in Ref. [l6| , Eq. ([8]) does not have a quartic term in A a ^ . 
This is because we neglected the color-chemical poten- 
tial /is and, in addition, we take into account only one 
dynamic gluonic field (see below). 

In the gluonic phase [l6[, the chromomagnetic insta- 
bility at A/<5/i < v2 triggers a non- vanishing vacuum 
expectation value of the spatial component of 



1 (A*-iA* 



(9) 



2 In the special case of a mixed phase, e.g., a color-electric field 
is generated around the boundary layer between the two phases 
and prevents the generation of a color-electric current across this 
layer. 



(For a simpler version of such a phenomenon, see also 
Ref. [25|.) Using the SO(3) rot rotational symmetry and 
the SU(2) C color symmetry, one can choose B = g{A & z ) ^ 
without loss of generality. Consequently, the non-zero 
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vacuum expectation value of B breaks S0(3) ro t, leaving 
only SO(2) rot [l^]. Furthermore, a non- vanishing B to- 
gether with A and /i e breaks the original symmetry of 
QCD down to 



U(lU®U(l) T 3®U(l) T s®SO(2) rc 



(10) 



where U(l) r 3^ is a subgroup of the SU(2)^//? chiral sym- 
metry and the charge Q is given by 



Q = Q f ®t c -tj 



y/3 



(11) 



with Qf — diag(|,— i) being the flavor matrix of the 
electric charges of quarks. 

The reduced symmetry of the ground state with B ^ 
allows for additional condensates, C = g{A\) j^O and 
D = g(Al) =/= 0. In fact, as discussed in Ref. fig , such 
condensates are required by the equation of motion. As 
discussed above, the gluonic field D is nothing but a color 
chemical potential fx^. The field C, on the other hand, in- 
duces electric superconductivity in the ground state and, 
therefore, is physically more interesting. However, in- 
cluding all three gluonic fields makes the analysis quite 
involved. In this work, we retain only the B field that is 
directly connected to the Meissner masses of gluons 4-7 
and, thus, is the most relevant field for the chromomag- 
netic instability. 

It is straightforward to show that the mass of the B 
field at B = (i.e., in the 2SC/g2SC phases) coincides 
with the Meissner screening masses of gluons of adjoint 
color 4-7 calculated in the hard-dense-loop (HDL) ap- 
proximation Q, 



Mi = 



d 2 v n 



dB 2 



J£_ 

6tt 2 



26fi 2 
A 2 



2 w5EAf 9(j;i _ A) 



(12) 



In order to derive this expression we neglected terms of 
order 0{fx 2 /A 2 ) and 0(A 2 /^ 2 ). 

In order to study the effect of the condensate field B 
on the free energy of the 2SC/g2SC phases, we calculate 
the difference of the thermodynamic potentials in a dense 
medium and in vacuum at the same value of B, 



VJA,B,Sfx,fx)-V g (0,B,0,0). 



(13) 



In a gauge theory, this subtraction in the one-loop free 
energy also takes care of the renormalization of the gauge 
coupling constant. As a result, the cutoff dependence of 
the free energy can be completely removed in this ap- 
proximation. 

Let us look at the free energy Q g in detail. Fig.[T]shows 
the free energy Q g (measured with respect to the normal 
phase at B = 0) as a function of A/Sfx. The results are 
plotted for fx = 500 MeV and 5fx = 80 MeV, with the 
diquark coupling chosen so that Aq = 132 MeV. Here 
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FIG. 1: The free energy (|13p measured with respect to the 
normal phase at B — 0, as a function of A/S/j, for B — MeV 
(solid), 50 MeV (dotted), 150 MeV (dashed), 250 MeV (dot- 
dashed), 350 MeV (short-dashed). The thin vertical line 
denotes the critical point of the chromomagnetic instability 
A/ Sfi = V2. We used fx = 500 MeV and 5li = 80 MeV, and 
the value of the coupling constant has been chosen so that 
Aq = 132 MeV. The contribution of electrons is not included. 
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FIG. 2: The free energy (|13[l measured with respect to the 
normal phase at B = 0, as a function of B in three regimes 
with different values of A/5fi. The solid lines correspond to 
the cases of A/S/j, — 0.8, 1.2, 1.6 (top to bottom). The dashed 
and dotted line denote the free energy at A/5fi = and at 
the critical point A/ 5 fx ~ v2, respectively. Other parameters 
are the same as in Fig. [1] 
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we do not restrict A/8fi to its physical value, determined 
by the stationary point of f2 ff (A), but treat it as a free 
parameter. 

Several important features of the free energy as a 
function of the B field arc evident from Fig. [TJ When 
A/S/i > v2) one can see that the free energy mono- 
tonically increases with B. (Strictly speaking, since our 
model reproduces the HDL result only up to terms of or- 
der 0(/i 2 /A 2 ) and 0(A 2 //z 2 ), the actual critical point is 
somewhat lower than V2.) In other words, Q g (B) has a 
global minimum at B — 0, and the 2SC phase is stable 
against gluon condensation in this regime. This is also 
clear from a different representation of the results, shown 
in Fig. H 

When A/ 8/i < a/2, on the other hand, we observe the 
onset of the chromomagnetic instability. For small B, the 
free energy first decreases with increasing B and then 
grows at larger B. This can be seen clearly in Fig. [2] 
The behavior of O s agrees well with Eq. (fT2"|) at small 
B. In this regime, the 2SC/g2SC phase is no longer the 
ground state. It is unstable with respect to the forma- 
tion of a non-zero B condensate, i.e., the so-called gluonic 
phase. The corresponding ground state is determined by 
the minimum of il g (B). The Meissner masses squared, 
which are given by the curvature of the free energy at the 
minimum, are non-negative in this state. [It is interest- 

I 



with the dimensionless parameter x\ being 
and 

q=\q], ?=-^4*>. (17) 

Note that the wave vector of the diquark condensate 
q is equivalent to a gauge field condensate (^4 8 ) in the 
case of single plane- wave LOFF pairing. In Eq. (fT5|) . the 
2SC/g2SC part of the free energy is obtained by taking 
the q — > limit. Also a non-zero color chemical potential 
fi$ has been neglected there. 

The free energy of a given phase can be computed by 
solving the gap equations, e.g., d^l/dA = and dtt/dq = 
0, and the neutrality condition dfl/dS^j, = 0. To simplify 
the calculations in the gluonic phase, we evaluate the 
free energy approximately as follows: (i) we obtain A* 
and Sfx* in the 2SC/g2SC phase by solving the coupled 



ing to note that, although the free energy in the normal 
phase, A/<5/i = 0, cf. dashed line in Fig. [3J increases 
with B, this increase is extremely slow over quite a large 
range of B values. The reason is that the quadratic term 
in the Taylor expansion of the free energy is vanishing, 
see Eq. (fT2")) .] Let us note that the results of Fig. [5] are 
in qualitative agreement with those shown in Fig. 2 of 
Ref. 0. 

From the results for the free energy it is clear that the 
gluonic phase resolves the chromomagnetic instability of 
the 2SC/g2SC phases. A neutral LOFF state is another 
candidate for the solution to the instability: it has been 
shown that such a state is free from the chromomagnetic 
instability, however, only in the weak-coupling regime 
(To| (for strong coupling, this is not the case \2l\). In or- 
der to determine the energetically most favored state, it 
is necessary to compare the free energies of the 2SC/g2SC 
phases, the neutral LOFF state and the gluonic phase. 

To this end we use the following approximation derived 
in Ref. [H| for the 2SC/g2SC phases and the neutral 
LOFF state: 

^ = ^2SC + ^ g 2SC/LOFF, (14) 

where the 2SC and g2SC/LOFF parts of the free energy 
are given by 



2/2 2 ) , (15a) 




(15b) 



I 

set of equations dfl/dA — and dQ/dSfi = 0, (ii) by 
using these solutions, we calculate n, g (B, A* ,8^,*) which 
is an approximate value for the free energy in the gluonic 
phase. For the densities of interests, B is at most of the 
order of 100 MeV, whereas q is of the order of tens of 
MeV. We performed a preliminary test of the quality of 
our approximation by varying B from to 300 MeV and 
computing the value of S[i necessary to ensure electric 
neutrality. We found that this value changes at most by 
10%. 

We illustrate the comparison of the free energies of 
all three phases in Fig. [3] (cf. Fig. 2 in Ref. [HI). We 
take \x = 400 MeV and choose the normal phase as a 
reference point for the free energy. In Ref. [2l|, it has 
been demonstrated that the neutral LOFF state is more 
stable than the 2SC/g2SC phases in the whole LOFF 
window 63 MeV < A < 137 MeV, which includes the 
entire g2SC window 92 MeV < A < 130 MeV. However, 
whereas the Meissner masses squared of gluons 4-7 in the 
weakly coupled neutral LOFF state are positive 0, [l(| , 
they remain negative in the intermediate and the strongly 
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FIG. 3: The free energy of the neutral 2SC/g2SC phase (solid 
line), the neutral LOFF state (dotted line), and the gluonic 
phase (dashed line) as a function of Ao. The three dots on 
the bold solid line (A = 92, 130, 177 MeV from left to 
right) denote the edge of the g2SC window with the nor- 
mal phase, the phase transition point between 2SC and g2SC 
phases (A = Sfi), and the critical point of the chromomag- 
netic instability (A = y/25)j,). The quark chemical potential 
is taken to be fi = 400 MeV. 



coupled regimes (at all values of Ao above 81 MeV) [2lj ] . 
In contrast, the gluonic phase removes the instability and 
is energetically favored over the 2SC/g2SC phases in the 
whole window in which the instability takes place (at all 
values of A/Sfi below y/2). The gluonic phase and the 
neutral LOFF state coexist in the region 92 MeV < Ao < 
137 MeV, but, as our results indicate, the gluonic phase 
is more stable than the neutral LOFF state in a wide 
region A > 103 MeV, which is close to the edge of the 
g2SC window with the normal phase. We argue therefore 
that the instability in the 2SC/g2SC phase is resolved by 
the formation of the gluonic phase. 

In summary, we explored the gluonic phase away from 
the critical point A/S/i — v2- We demonstrated that 
the energetically favored state of a neutral two-flavor 
color superconductor is not the 2SC/g2SC phase in the 
intermediate- and strong-coupling regimes but the glu- 
onic phase in which the dynamic gluonic field B = g(A®} 
acquires a vacuum expectation value. In particular, the 
whole g2SC phase is replaced by the gluonic phase which 
is chromomagnetically stable. 



We also compared the free energies of the 2SC/g2SC 
phase, the neutral LOFF state, and the gluonic phase. 
We found that the gluonic phase is energetically favored 
in the intermediate- and strong-coupling regimes. The 
encouraging results of this analysis should be further im- 
proved in the future by (i) taking into account the most 
general ansatz for the gauge-field configuration in the glu- 
onic phase, (ii) by calculating the free energy in a self- 
consistent manner. We already performed a preliminary 
investigation including the effect of the gluon field D, 
responsible for enforcing Gauss' law, and found that, for 
Ao < 130 MeV, the additional cost in the free energy is of 
order 0.01 MeV/fm 3 , and thus negligible. For Ao > 130 
MeV, however, the effect of the D field could be an order 
of magnitude larger. 

It is appropriate to mention that the instability re- 
lated to the 8th gluon was not studied in the present 
work. In this sense, the neutral LOFF state is appealing, 
because the Meissner mass of the 8th gluon is automati- 
cally zero in this state. It should also be mentioned that 
in the strongly coupled LOFF state 0, [ll| the longitudi- 
nal Meissner mass squared of the 8th gluon is negative. 
Although this instability was not addressed in this work, 
it is unlikely, however, that the LOFF state is energeti- 
cally favored in the strong-coupling regime. 



Note added in proof 

It has recently been demonstrated that, in the three- 
flavor case, realistic crystal structures are more robust 
than a single plane-wave LOFF state 
flavor case, Bowers and Rajagopal 



22j. In the two- 



already indicated 
that a LOFF state with multiple plane waves would have 
a lower free energy than that with a single plane wave. 
The result shown in Fig. [3]would be altered by the inclu- 
sion of crystal structures with more plane waves. This is 
therefore an important project that needs to be addressed 
in future work. 
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